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13.  ABSTRACT 


integrals  with  integrands  of  the  form  )  £(t)  are  considered  for  A  -  -  ' 

ll(t)  oscIiiat°ry  for  large  argument.  It  is  shown  that  the  set  of  critical 

po'.uis  for  such  integrals  includes  zeros  of  the  phase  function  *  as  well  as  all 

of  those  that  arise  in  the  analysis  of  the  standard  integrals  of  Fourier  type; 

for  lIi.  rial  case  '-here  H(t)  -  oxp(it)  .  The  contribution  to  the 

asymptotic  expansion  from  each  type  of  critical  point  is  derived.  In  particular, 

a  formula  is  obtained  which  rener  1 1  i  fU.,  of  .i  •  , 

s'  ncrauz.  s  the  stationary  phase  lormula  associated 

with  Fourier  type  integrals. 


D D  ,1473 


1.  Introduction. 


In  the  method  of  stationary  phase,  one  is  concerned 
with  the  asymptotic  expansion,  as  X  -  •,  of  functions 
defined  by  integrals  of  the  form 

,b 

(1.1)  I(X)  -  I  exp  {iXcr>(t))f(t)dt. 

J  a 

We  shall  assume  that  the  details  of  this  method  are 
familiar  to  the  reader  and  need  not  be  discussed. 

Our  concern  here  shall  be  with  the  asymptotic 
expansion  of  integrals  of  the  form 

„b 

(1.2)  I(X)  =  f  H(  Xcr(t)  )f  (t)dt 

Ja 

in  the  case  where  the  kernel  function  H(t)  is  oscillatory 

for  both  large  positive  and  large  negative  arguments. 

More  precisely,  we  assume  that,  as  t  -  +*», 
r  Kv/ 6  .  <r  N(m)  _r 

(1.3)  H(t)  ~expji  £  V*  /X  ^  c^t  “(log  t)n. 

1=0  m=0  n=0 

Here  each  b^  Is  real,  N(m)  is  finite  for  all  m  and 
(Re(rm)]  is  a  monotonically  increasing  sequence  with 
limit  +«>.  In  the  limit  t  -  -®,  we  assume  that  an 
expansion  of  the  form  (1.3)  holds  with  t  replaced  by 
|t|  and  with,  in  general,  different  constants. 

Clearly  the  Fourier  kernel  exp{it!  is  a  special 
case  of  the  general  kernel  we  propose  to  study.  Thus 
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we  should  expect  to  recover  from  our  asymptotic  analysis 
of  (1.2)  the  stationary  phase  results  valid  for  (l.l). 

We  should  point  out  that  there  are  functions,  such 
as  the  Airy  function  Ai(t),  which  are  oscillatory  in 
one  of  the  limits  t  -  ±»  and  exponential  in  the  other. 
It  will  be  apparent  that  integrals  (1.2)  with  sue;  iunc- 
tions  as  kernels  can  also  be  treated  by  the  methods  to 
be  developed  below. 

The  method  itself  involves  applications  and 
generalizations  of  an  asymptotic  technique  recently 
developed  by  Handelsman  and  Lew  [1],  [2],  [3].  This 
technique  makes  heavy  use  of  the  Mellin  transform  whose 
relevant  properties  are  discussed  below.  In  references 
[  1  ],  [  2  ]  and  [  3  ]  only  the  case  cp(t)  =  t  is  treated. 
Here,  however,  we  shall  consider  more  general  cp  and, 
in  particular,  shall  allow  «o  to  be  non-monotonic. 

As  is  well-known,  to  derive  the  asymptotic  expansion 
of  (1.1)  one  must  first  identify  its  set  of  critical 
points.  These  include  end  points  of  integration, 
stationary  points  of  ra  and  points  where  either  cr>  or 
g  feils  to  be  infinitely  differentiable.  We  shall  show 
that,  in  addition  to  all  of  the  above,  the  set  of 
critical  points  for  (1.2)  also  includes  the  zeros  of 
Indeed,  this  is  one  of  the  main  results  of  this  paper. 

To  gain  some  insight  into  the  critical  nature  of 
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the  zeros  of  ©,  suppose  that  (1.3)  is  only  an  asymptotic 
result  and,  in  particular,  is  not  valid  near  t  *  0. 

Then,  no  matter  how  large  X  is,  there  always  exists  a 
neighborhood  of  each  zero  of  ©  throughout  which  H( X©) 
is  not  asymptotically  described  by  our  assumed  asymptotic 
forms.  In  other  words,  the  asymptotic  expansion  of  H(X©), 
as  X  -  undergoes  a  drastic  change  as  t  passes 
through  any  of  these  neighborhoods.  For  this  reason  we 
can  think  of  these  neighborhoods  as  "boundary  layer 
regions."  It  is  certainly  reasonable  that  the  rapid 
change  in  the  asymptotic  behavior  of  H( X©)  as  t 
passes  through  a  boundary  layer  region  will  affect  the 
asymptotic  expansion  of  I. 

In  the  light  of  the  above  argument  we  can  under¬ 
stand  why  the  zeros  of  ©  are  not  critical  pointB  for 
(1.1).  Indeed,  the  Fourier  kernel  H(t)  =  exp(it), 
has  an  asymptotic  expansion  as  t  -  +®  of  the  form 
(1.3).  This  expansion,  however,  holds  for  all  t  so 
that  there  are  no  boundary  layer  regions  of  the  type 
just  described. 

In  the  following  section  we  reduce  our  problem  to 
the  study  of  certain  integrals  of  canonical  type.  In 
Section  3,  we  consider  some  results  concerning  Mellin 
transforms  that  are  needed  to  implement  our  methods. 

Finally  the  desired  asymptotic  expansion  of  I  is. 
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obtained  in  Sections  4  and  5. 

2.  Reduction  to  Canonical  Integrals. 

Because  there  are  many  possible  critical  points 
for  integrals  of  the  form  (1.2)  with  H  an  oscillatory 
kernel,  it  is  convenient  to  have  a  means  for  isolating 
them  so  that  their  contributions  to  the  asymptotic 
expansion  of  I  can  be  studied  separately.  This  can 
be  accomplished  by  using  neutralizer  functions.  These 
were  first  introduced  by  Van  der  Corput  [4  3,  and  we 
shall  assume  that  their  basic  properties  are  familiar 
to  the  reader.  The  net  effect  of  the  neutralization 
process  is  to  reduce  the  asymptotic  analysis  of  (1.5) 
to  the  study  of  a  sum  of  Integrals  each  having  exactly 
one  critical  point  either  as  an  upper  or  as  a  lower 
end  point  of  integration. 

Suppose  first  that  t  =  tQ  is  a  critical  point 
at  which  ©  is  non-zero.  After  neutralization,  t  =  tQ 
will  appear  as  either  an  upper  or  lower  end  point  of 
integration  in  at  most  two  of  the  integrals  to  be 
asymptotically  evaluated.  To  obtain  the  corresponding 
contributions  to  the  asymptotic  expansion  of  I,  one 
need  only  replace  H(Vi)  by  the  appropriate  asymptotic 
expansion  and  integrate  the  resulting  series  term  by 
term.  Thus,  finding  these  contributions  is  reduced  to 
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the  asymptotic  evaluation  of  many  integrals  of  the  form 
(1.1).  We  have,  therefore,  that  the  only  critical  points 
which  require  non-standard  methods  of  analysis  are  the 
zeros  of  o. 

As  a  result  of  the  above  discussion,  we  shall  focus 
our  attention  on  obtaining  the  contribution  to  the 
asymptotic  expansion  of  (1.2)  corresponding  to  a  given 
zero  of  q.  If  we  denote  this  zero  by  t  =  c  and  the 
contribution  by  IC(X),  then  after  neutralization  we 
find  that 


(2.1) 

where 

(2.2) 


!<.(*>  =  Ic+(>l)  +  Jc  (A) 


Ic  (*)  =  J  H(Xcc)gc  (t)dt 


and 


(2.?) 


-D 

!c  (X)  =  j  H(Mgc  (t)dt. 


Here  gc  (t)  vanishes  for  t  <  a  <  c  with  a  chosen 

so  that  ©  vanishes  in  fa,c]  only  at  t  =  c  and  ©• 

does  not  vanish  in  [a,c).  Furthermore,  gc  (t)  =  g 

in  some  small  half  neighborhood  of  t  =  c  . 

Similarly,  g  (t)  vanishes  for  c  <  3  <  t  with 
G+ 

3  chosen  so  that  ©  vanishes  in  [c,p]  only  at  t  =  c 


and  ©'  does  not  vanish  in  (c,3).  Finally,  g  =  g 

c+ 
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in  some  small  half  neighborhood  of  t  =  c+.  Of  course 
if  c  coincides  with  one  of  the  end  points  of  inte¬ 
gration  in  (1.2),  then  only  one  of  the  integrals  1  , 

I  is  non-zero. 
c+ 

Suppose  now  that  as  t  -  c_,  ep  has  an  asymptotic 
expansion  whose  leading  term  is  given  by 

(2.4)  *  ~  Y0(c-t)V°,  v0  >  0. 

If  in  (2.2)  we  introduce  the  new  variable  of  integration 

(2.5)  s  -  uc  ®(t),  uc  =  sgn  yQ 
then  we  can  write 

(*<** 

(2.6)  I  (A)  «  )  H(Xm  s)G  (s)ds. 

c-  *■  0  - 

Here 

(2.7)  G_  =  gc  ('t(s))§|- 

which  we  note  vanishes  for  s  >  uc  ®(0). 

Similarly,  if  we  assume  that  as  t  -  c+ 

(2.8)  n(t)  ~  ^Q(t-c)P°,  o0  >  0, 
and  set 

(2.9)  £  =  u_  ©,  u  =  sgn  vn, 

S  c+  u 

then  I  can  be  written 
c 
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(2.10)  Ic  (X)  =  I  H(  Xp  s)G  (s)ds. 

t'+  J  0  c+  + 

Here 

(2.11)  G+(s)  =  gc  (t(s))§| 

which  vanishes  for  s  >  u  o(a). 

c+ 

Thus  we  have  reduced  our  problem  to  the  study  of 
the  two  canonical  integrals 

r00 

(2.12)  I  ( X)  =  |  H(±Xs )G(s )ds 

*  “0 

where  the  kernel  H(t)  is  oscillatory  in  each  of  the 
limits  t  -  *«  and  G(s)  vanishes  for  s  outside  of 
some  finite  interval. 

3.  Results  on  Mellln  Transforms. 

As  we  shall  see,  Mellin  transforms  play  an  important 
role  in  our  asymptotic  development.  Indeed,  one  might 
anticipate  this  upon  observing  that  each  of  the  canonical 
integrals  (2.12)  can  be  expressed  as  a  Mellin  convolution 

[3  ]. 

The  Mellin  transform  of  a  function  f(s)  is  defined 

by 

09 

(3.1)  M[f;z]  =  f  f(s)sz-1ds,  z  «  x+iy 

J0 

when  this  integral  exists.  Furthermore  if  f  is  such  that 
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(3.2)  f(s)  =  0(sp),  s  -  0+ 

f  (  S  )  =  0(  S  )  ,  S  — *  +«• « 

Then  M[f;z]  converges  and  is  holomorphic  in  the  strip 
(3-3)  -p  <  Re(z)  =  x  <  r. 

Also,  within  this  strip  lim  |M[f;x+iy]|  «=  0. 

I  y  I  — 

If  in  (3.2)  -p  >  r,  then  M[f;z]  does  not  exist 
in  the  ordinary  sense.  Nevertheless  with  certain  addi¬ 
tional  assumptions  on  f  one  can  define  M[f;z]  in  a 
generalized  sense  (  3  ] .  We  shall  not  need  to  appeal  to 
this  extension  in  this  paper  however. 

From  our  point  of  view,  there  are  two  results  con¬ 
cerning  Mellin  transforms  that  are  of  special  significance. 
The  first  is  the  simple  relation 

(3.-0  M[f(As);z]  =  VZf  f  (s  )sZ_1ds  =  A_ZM(f;z]. 

'  0 

The  second  involves  integrals  of  the  form 

(3-5)  J  =  r  f(s)h(s)ds. 

'  0 

Indeed,  suppose  that  M[f;z]  and  M[h;z]  are  holo¬ 
morphic  in  overlapping  vertical  strips.  (This  will 
always  be  the  case  if  J  is  absolutely  convergent.) 

If  Re(z)  =•  c  lies  in  the  common  strip  of  analyticity, 

tr- n  we  na  e 
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r®  i  (*c+i® 

(3.6)  f (s)h(s)ds  =  *^-1  M[h;z]M[f;l-z)dz 
J0 

which  is  Parseval's  Theorem  for  Mellin  transforms. 

Upon  combining  these  last  two  results,  we  find 
that  our  canonical  integrals  (2.12)  have  the 
rep  re  s  ent  at ions 

,  c  +i® 

(3.7)  1  M  "  2nT!  X“zM[H(±s);z]M[G(s);l-z]dz. 

c^-i® 

Here  Re(z)  *  c  lies  in  the  strip  of  analyticity  of 
the  integrand.  We  note  that  the  total  dependence  of 
the  integrand  on  X  is  contained  in  the  factor  >TZ . 

Our  plan  is  to  push  the  contour  of  integration  in 

(3.7)  to  the  right,  apply  Cauchy's  integral  theorem, 
and  derive  thereby  an  asymptotic  expansion  of  1*(X) 
as  a  residue  series.  In  order  to  accomplish  this  we 
must  obtain  certain  information  about  the  analytic 
continuations  of  the  functiona  M[H(±s);z]  and 
M[G(s);l-z]  into  the  right  half  plane.  Specifically, 
we  must  locate  and  classify  the  singularities  of  these 
continuations,  and  we  must  estimate  their  behavior  as 
z  -  ®  along  vertical  lines. 

For  oscillatory  functions,  the  required  information 
is  contained  in 

Lemma  1.  Suppose  that  H(s)  is  locally  integrable  on 
(0,®)  and  satisfies  (3.2)  with  -p  <  r.  Suppose  further 
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that,  as  s  -•  -h»,  H(s)  has  an  asymptotic  expansion  of 
the  form  (1.3)  In  which  event  -r  *  rQ.  Then  MfHjz) 
can  be  continued  Into  the  right  half  plane  Re(z)  >  -p 
as  a  holomorphic  function.  Furthermore,  in  this  right 
half  plane 

/  (x-Re(rn))/v-i\ 

(3.8)  | M[H;z] j  =  o()y|  0 

as  |  y |  -  ®. 

Proof:  The  proof  of  this  lemma  is  given  in  the  Appendix. 

As  an  example  and  to  illustrate  the  sharpness  of 
the  estimate  (3.8),  let  us  consider  the  function 
H  =  exp{isv)  whose  Mellin  transform  is  given  by  [ 5  ] 

(3-9)  M[exp(isv);z]  =  exp^^^f^),  0  <  x  <  1. 

In  this  case  rQ  =  0,  and  the  analytic  continuation 
into  Re(z)  >  1  is  explicit.  We  note  that  r(^  is 
analytic  in  Re(z)  >  0.  The  estimate  (3-8)  follows 
from  the  known  asymptotic  expansion  of  the  gamma 
function 

*4 

(3.10)  r(z)  ~  j y |  exp[-TTjy|/2],  |y|  -  «. 

As  we  have  indicated  above,  our  plan  is  to  push 
the  contour  of  integration  in  (3-7)  to  the  right.  In 
order  to  accomplish  this  we  must  of  course  determine 
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the  analytic  continuation  of  M[G;l-zJ  into  a  right 
half  plane.  Let  us  assume  for  the  present  that  this 
has  been  done  and  M[G;l-zJ  is  a  meromorphic  function. 

Then  to  Justify  the  displacement  of  the  contour  to  the 
line  Re(z)  «=  k  >  c±,  one  must  still  show  that 

(3.11)  lim  M[H(±s);z]M[G(s);l-z]  =  0,  c±  <  x  <  k. 

I  y  I 

The  estimate  (3-8)  implies  an  algebraic  growth  of 
M[H(±s);z]  in  this  limit  which  worsens  with  increasing 
x.  This  growth  must  therefore  be  compensated  by  a 
commensurate  decay  of  the  analytic  continuation  of 
M[G(s);l-z].  In  the  following  sequence  of  lemmas,  we 
shall  establish  sufficient  conditions  for  such  decay. 

Lemma  2.  Let  G(s)  be  q  times  continuously  differentiable 
on  (0,®).  Let  G^q+1^(s)  be  piecewise  continuous  on 
[0,k]  and  continuous  for  s  >  k.  Finally  suppose  that 
there  exists  a  real  number  Xq  such  that  for  all  x  >  Xq, 
(sSsf  (sXG(s))  vanishes,  as  s  -  0+,  for  p  =  0,l,2,...,q 
and,  as  s  -  for  p  =  0,1,2, ...  ,q+l.  Then,  as  ly|  -»  ®, 

(5-12)  M[G;z]  ®  0(|yrq'1) 

for  all  x  >  xQ. 

Proof :  The  proof  of  this  lemma  is  given  in  the  Appendix. 
Remarks :  The  hypotheses  of  Lemma  2  simply  provide 
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sufficient  information  to  allow  for  the  estimation  of 
M[G  jz]  via  integration  by  parts.  Furthermore  the 
assumptions  on  G  imply  that  M[G;z]  is  holomorphic 
in  — Xq  <  Re(z). 

The  next  two  lemmas  follow  from  analogous  results 
for  Fourier  transforms.  Their  proofs  will  be  omitted 
here,  but  can  be  constructed  from  the  corresponding 
proofs  in  Titchmarsh  { 6  ] ■ 

Lemma  3.  Let  G(s)  satisfy  the  conditions  of  Lemma  2 
except  now  replace  the  condition  on  G^q+1^(s)  by  the 
assumption  that  (sA)(sXg(s))  is  of  bounded  total 
variation.  Then 

(3.13)  MfGjz]  »  0(|y|"^-1), 
as  |  y |  -  <*>,  for  all  x  >  xQ. 

Lemma  4.  Let  G(s)  satisfy  the  conditions  of  Lemma  2 
except  now  replace  the  condition  on  G^+^(s)  by  the 
assumption  that  is  Hftlder  continuous  of 

order  y  on  [0,k]  and  of  bounded  total  variation  for 
s  >  k.  Then 

(3.14)  M[G;z]  =  0(|yrq"7) 
as  jy|  -  »•,  for  all  x  >  xQ. 
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Lemmas  2-4  yield  estimates  on  the  decay  of  M[G;z] 
in  its  region  of  absolute  convergence  or  equivalently 
in  its  region  of  analyticity.  We  now  wish  to  obtain 
analogous  information  outside  of  this  region.  Ao  we 
shall  soon  see,  the  analytic  continuation  of  M[G;z] 
to  the  left,  (and  hence  of  M[G;l-z]  to  the  right), 
depends  to  a  large  extent  on  the  nature  of  G(s)  near 
s  *  0+.  Indeed  we  have  the  following  result  due  to 
Handelsman  and  Lew  [ 3  ] : 


Lemma  5.  Suppose  that  M[G;z]  is  holomorphic  in  the 
region  -a  <  Re(z)  <  S,  and  that,  as  s  -  0+, 


"S? 


V8"”'10* 8  >" 


(3.15)  O(s)  ~  £ 

m=0  n-0 

with  Re(affl)  t  ®  and  N(m)  finite  for  each  m.  Then, 

a  -  Re(aQ)  in  (3.2)  and  M[G;z]  can  be  continued  into 

Re(z)  <_  -  Re(aQ)  as  a  meroinorphic  function  with  poles 

at  the  points  z  =  -am.  Moreover,  about  these  points, 

M[G;z)  has  a  Laurent  expansion  with  singular  part 

N(m) 

.  y  d  -  x^i).. 

L  mn#  _  \n+i 
n=0  (  2  V 


Remark.  We  note  that  when  N(rrO  =  0  for  each  m,  i.e., 
when  no  logaritnms  appear  in  the  expansion  (3.15)>  all 
of  the  poles  in  the  analytic  continuation  of  M[G  ;z]  to 
the  left  are  simple. 
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If  we  combine  the  results  of  Lemmas  1  and  5,  then 
we  can  conclude  that  when  G(s)  has  an  expansion,  as 
s  -  0+,  of  the  form  (3.15),  ell  of  the  singularities 
of  the  analytic  continuation  of  M[H(±s);z],  M[G;l-z] 
into  the  right  half  plane  are  determined  by  the 
exponents  am«  Moreover,  these  singularities  are  poles 
so  that  our  proposed  deformation  of  contour  will  indeed 
yield  a  residue  series  for  I*(*).  We  must  still  esti¬ 
mate  M[G;l-z]  as  J y J  -  «  in  order  to  Justify  the 
deformation.  For  this  purpose  we  now  state 

Lemma  6.  Let  G(s)  satisfy  the  smoothness  conditions  of 
Lemma  2.  Also,  let  vanish,  as  s  -  +», 

for  p  =  0,1,..., q+1  and  x  >  l-Re(aQ)  «=  xQ.  Finally 
suppose  that  (3.15)  holds  and  that  the  asymptotic 
expansion  of  dm^(s),  m  ■  0,...,q+l,  as  s  -  0+, 
is  obtained  by  successively  differentiating  (3.15  ) 
term  by  term.  Then 

(3.16)  M[G  ; z ]  -  0(  |y|  "<*"1) ,  |y| 

for  all  x.  Here  by  M[G;z]  we  mean  the  analytic  con¬ 
tinuation  of  this  Mellin  transform  into  the  entire  z-plane. 

Proof:  The  proof  of  this  lemma  is  given  in  the  Appendix. 


Corollary .  If,  in  Lemma  6,  the  stated  conditions  hold 
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~x* 

for  all  q,  then  M[G;z]  *  0(|y|  )  for  all  r  and 

all  x. 

We  remark  that,  if  In  Lemma  6,  the  smoothness 
conditions  of  Lemma  2  are  replaced  by  those  of  either 
Lemma  3  or  Lemma  4,  then  the  corresponding  changes 
must  be  made  in  the  estimate  (3.16).  Nevertheless,  one 
still  finds  that  the  results  obtained  are  valid  for  the 
analytic  continuation  of  the  Mellin  transform  into  the 
entire  complex  plane. 

To  illustrate  some  of  the  results  obtained  above 
let  us  consider  an  explicit  example.  Indeed,  suppose 
that 

(3.17)  G(s)  -  sae"s 

which  satisfies  the  conditions  of  Lemma  6  with  aQ  =  a 
and  q  *  o'.  Then  the  corollary  predicts  that,  as  |y|  -  ®, 
M[G;z)  decays  faster  than  any  power  of  |y|.  For  this 
example  we  have  the  explicit  result 

M[G;z]  -  f(a+z)  =  ofexpf^El)),  |y|  -  ~ 

which  agrees  with  this  prediction.  Furthermore  we  have 
from  known  properties  of  the  gamma  function  that  the 
analytic  continuation  of  r(a+z)  into  Re(z)  <  -a  has 
simple  poles  at  the  points  z  =  -(a+m),  m  «=  0,1,2,..., 
with  corresponding  singular  parts 
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(3.18)  ( -l)m/m! ( z+a+m) . 

As  is  readily  seen,  this  last  result  is  in  agreement 
with  that  predicted  by  Lemma  5. 


4.  Asymptotic  Expansion  of  I±(X). 

By  using  the  theory  of  Mellin  transforms  developed 
in  the  previous  section,  we  shall  now  derive  asymptotic 
expansions  for  the  two  canonical  integrals  (2.11).  We 
first  note  that,  if  H(±s)  has  the  asymptotic  expansion 
(1.3)  as  s  -  and  G(s)  has  the  asymptotic  expansion 
(3.13)  as  s  -  0+,  then  in  (3.7) 


(4.1) 


c±  <  Min(Re(rQ),Re(l+a0)) 


since  this  is  the  right  limit  of  the  common  strip  of 
analytic ity  of  the  integrand  in  that  equation. 

We  state  the  main  result  concerning  the  asymptotic 
behavior  of  I*(X)  in 

Theorem  1.  Let  G(s)  satisfy  the  conditions  of  Lemma  6 
and  H(±s)  satisfy  the  conditions  of  Lemma  1.  Then 

(am+1) 


N(m) 


r 

L 


,  d  X 

l,  mn 


(4.2)  I*(X)  = 

Re(am+l)<k  n=0 

n  /  \  c 

X(;)(-log  X)JM(n-j)[H(is);l+am]  +  (f(X;k). 

j=0 


Here 
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(4.3) 


,  „k+i*  _ 

X  zM[H(±s ) ;z]M[G(s ) ;l-z) 

k*i# 

o(X  k),  X  —  «D 


and 

(4.4)  k  <  v(Q+f)  +  Re(rQ) 
where  k  /  Re(am)+1  for  any  m. 

Proof :  In  the  exact  representation  (5.7)  we  displace 
the  vertical  contour  of  integration  to  the  line 
Re(z)  =  k  >  c±.  We  note  that  by  Lemmas  1  and  5  the 
analytic  continuation  of  the  integrand  in  (3-7)  into 
the  right  half  plane  is  a  meromorphic  function  with 

poles  at  the  points  z  =  am+l,  m  =  0,1 .  Indeed, 

we  find,  upon  formally  applying  Cauchy's  integral 
theorem,  that  (4.2)  is  valid.  Thus  to  complete  the 
proof  of  the  theorem  we  need  only  justify  the  displace¬ 
ment  Itself  and  establish  the  error  estimate  given  by 
(4.3)  and  (4.4). 

It  follows  from  (3.8)  and  (3.16)  that 

(4.5)  M[H(±s);x+iy]M[G(s);l-x-iy]  =  0(|yfe(x)),  |y|  - 

Here 

e(x)  =  q  +  \  -  (x-Re(rQ))/v. 


(4.6) 
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Thus  we  can  displace  the  contour  to  the  line  Re(z)  *  k 
so  long  as  e(k)  >  0,  i.e.,  so  long  as 

(4.7)  k  <  v(q  +  |-)  +  Re(rQ). 

With  (4.7)  satisfied  we  have  that  (^(X;k)  exists. 

The  estimate  (4.5)  need  not  hold,  however.  We  note  that 
(5(*;k)  can  be  viewed  as  a  Fourier  transform  with 
respect  to  log  X.  Indeed,  we  have 

(4.8)  Xk£(X;k) 

=  exp{-iy  log  X]M[H(±s )  ;k+iy]M[G(s )  ;l-k-iy)dy. 

—  OB 

Suppose  now  that  (4.4)  is  satisfied  which,  in  turn, 
implies  that  e(k)  >  1.  Then  by  applying  the  Riemann 
Lebesgue  lemma,  we  find  that,  as  X  -  <*,  the  right  side 
of  (4.8)  is  o(l),  and  the  estimate  (4.5)  follows. 


Corollary.  Let  H(+s)  satisfy  the  hypotheses  of  Lemma  1 

and  let  G(a)  satisfy  the  hypotheses  of  Lemma  6  with 

q*=  ».  Then  the  infinite  expansion 
<»  N(m)  _/_  \  n  ,  ^ 

/i.  -v  m  V  /  r t _ 5\Jw(n“.l)r„/ 


<’!-9 >J*(»)  “  I  I>n,n 
m=0  n=0 


\  (j)(-log  X)JM(n"j)[H(±s);l+am] 


holds  as  X  -  ®. 


Proof :  It  follows  from  the  corollary  to  Lemma  6,  that 


in  this  case 
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(4.10)  M[H(±a) jx+ly)M[G(s) jl-x-iy]  =o(|y|"r),  (y|  - 

for  all  r  and  all  x.  Hence  we  can  let  k  go  to  +® 
in  (4.2)  and  (4.3)  to  obtain  the  desired  result. 

We  note  that  when,  in  (3.15),  d^  «  0  for  n  >  1 
and  all  m,  i.e.,  when  no  logarithms  appear  in  the 
asymptotic  expansion  of  G(s),  as  s  -  0+,  the  asymptotic 
expansion  (4.9)  reduces  to 

“  -(a  +1) 

(4.11)  I(X)  ~  X  dmOX  M[H(±s);l+am]. 

m*=0 

In  the  proof  of  Theorem  1,  the  Riemann  Lebesgue 
lemma  was  applied  to  estimate  Xk£(X;k)  as  X  -  ®. 
Recently,  Bleistein,  Handelsman  and  Lew  [7],  have 
obtained  a  generalization  of  this  lemma  which,  under 
slightly  more  restrictive  assumptions  can  be  used  to 
improve  the  error  estimate  found  in  Theorem  1.  Indeed, 
we  have 

Theorem  2.  Let  the  hypotheses  of  Theorem  1  be  satisfied 
so  that  (4.5)  holds  with  e(x)  defined  by  (4.6).  Suppose 
further  that 

(4.12)  M[H(±s )  ;x+iy]M[ G(s ) ;  1-x-iy]  ~  c0exp{iayu]|yre^x^, 

as  ( y !  -  where  a  and  u  are  any  real  numbers. 

Then  (4.2)  and  (4.3)  hold  with 
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(4.13)  k  <  v(q  +1 )  +  Re(rQ),  k  /  Re(am+l),  m  =  0,1,2,... 

Proof ;  It  follows  from  Theorem  1  that  all  we  need  show 
Is  that  (4.3)  holds  for  all  k  satisfying  (4.13).  Thus 
consider  Xk£(X;k)  as  given  by  (4.8).  The  results  of 
reference  [7  ]  show  that  whenever  e(k)  >  0, 

(4.14)  limXk£(X;k)  =  0, 

X-«° 

except  possibly  when  1  <  u  <  2  in  (4.12).  It  is 
further  shown  that  (4.14)  still  holds  with  1  <  p  <  2, 
so  long  as 

(4.15)  0  <  e(k)-(l-£)  =  q+l-(k-Re(r0))/v  +  |-(l-£). 

Since  1  3  p,  it  is  clear  (4.15)  is  satisfied  whenever 
(4.13)  holds.  This  completes  the  proof. 

We  wish  to  emphasize  that  the  hypotheses  of  Theorem  2 
differ  from  those  of  Theorem  1  only  in  that  the  former 
includes  an  additional  assumption  concerning  the  oscil¬ 
latory  behavior  of  M[ H( ±s ) j z ]M[ G(s ) ;l-z]  in  the  limit 
! y I  -  o'.  It  is  this  more  specific  information  that  allows 
us  to  apply  the  results  of  reference  [7  ]  and  thereby 
extend  the  validity  of  (4.3)  to  the  region  (4.4). 

Our  concern,  of  course,  is  ultimately  with  the 
integrals  I  ( x)  from  which  the  canonical  integrals 
I  ( X)  were  directly  derived.  We  recall  that  s  =  0  in 
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I*(X)  corresponds  to  t  =  c  in  I  (X)  where  c  is 

c± 

a  point  in  the  original  domain  of  integration  at  which 
the  phase  function  ®  vanishes,  if  we  assume  that  the 
conditions  of  the  corollary  to  Theorem  1  hold,  then  we 
find  that  s  =  0  is  the  only  critical  point  for  i+(X). 
Thus  the  infinite  expansion  (4.9)  can  be  used  to  obtain 
the  contribution  to  the  asymptotic  expansion  of  (1.2) 
corresponding  to  a  zero  of  ??.  We  must  point  out,  how¬ 
ever,  that  we  have  not,  as  yet,  established  the  critical 
nature  of  t  =  c.  This  can  be  done  by  explicitly 

obtaining  the  expansions  of  I  (x)  and  adding  them. 

c± 

Only  if  the  resulting  sum  is  non-trivial  can  we  conclude 
that  t  =  c  is  a  critical  point  for  I(X).  We  shall 
investigate  this  point  further  in  the  following  section 
along  with  some  illustrative  examples. 

5*  Explicit  Results  and  Examples. 

We  wish  now  to  determine  the  contribution  to  the 
asymptotic  expansion  of  (1.2)  corresponding  to  an 
interior  zero  of  Moreover,  we  want  to  express  this 

result  explicitly  in  terms  of  the  original  functions  r* 
and  g.  In  principle,  we  could,  by  using  the  results 
of  the  previous  section,  find  as  many  terms  of  this 
contribution  as  desired.  The  computations,  however, 
become  exceedingly  awkward  as  the  number  of  terms  increases 


22 


and  hence,  for  the  most  part,  we  shall  be  content  here 
with  obtaining  expansions  to  leading  order  only.  We 
shall  assume  throughout  this  section  that  the  functions 
«r>  and  g  are  sufficiently  smooth  so  that  either 
Theorem  1  or  2  can  be  applied  to  obtain  the  expansions 
to  the  orders  stated. 

Let  us  suppose  that  in  (1.2)  cr(c)  =  0  with 
a  <  c  <  b.  If,  as  in  Section  2,  we  denote  the  contri¬ 
bution  corresponding  to  t  =  c  by  Ic(X),  then  we  have 

(5.1)  Ic  -  I  (X)  +  Ic  (X) 

with  I  (x)  defined  by  (2.10)  and  (2.6)  respectively. 
We  now  assume  that,  as  t  -  c+, 

(5.2)  g(t)  ~  g+(t-c)  +  ,  w(t)  ~  T>0(t-c)  0 

0o-l 

C'(t)  ~  P0"0(t-c)  ,  o0  >  0. 


as  s 


0+.  Hence  we  find  from  (4.11)  that,  in  this  case 
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(5-5)  I.  (X) 
c+ 


as  \  -  ®.  Here  u+  =  sgn  nQ. 
Similarly,  if  as  t  -  c 


(5-6) 


u>  -1 

g(t)  ~  g_(c-t)  "  , 

cp' - vQyQ(c-t) 


w 

v0'1 


70(ct)V°, 

vQ  >  0. 


Then  we  find 

g  /  i  \o>  /vn  r  o>  . 

<5-7>  V11)  *  VW)  MrH(u-S,;^ 

Here  \i_  =  sgn  7Q. 

Upon  adding  (5*5)  and  (5*7)  we  obtain  the  desired 
contribution  from  t  =  c  to  leading  order.  In  most 
instances  the  constants  in  the  assumed  expansions  (5.2) 
are  closely  related  to  the  corresponding  constants  in 
(5.6).  Two  cases  are  worthy  of  special  consideration. 
Suppose  first  that  g  is  continuous  and  non-zero  at 
t  =  c  so  that 


♦  oCa--"0). 


(5-8)  g+  =  g_  =  g(c),  a>+  =  =  1. 

Suppose  further  that  cp  is  differentiable  at  t  =  c 
with  fp'  (c)  /  0.  Then 

(5-9)  n0  =  -70  =  <p' (c),  oQ  =  Vq  =  1. 

It  follows  from  (5. 5 -5. 9)  that 
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(5.10)  IC(X)  ~  |<pf||]|X(M[H(s);l]  +  M[H(-s)  ;1] } . 

Let  us  now  suppose  that  relations  (5.8)  hold,  but 
that  c?  has  a  simple  stationary  point  at  t  =  c.  Then 


(5.11)  n0  -  70  =  cp"  ( c  )/2 ,  PQ  =  vQ  =  1 


so  that,  in  this  case,  we  have 

(5.12)  Ic(A)  ~  g(  c  )  fy]  ^)1/2M[h(  sgn  cp"(c)s);|"|. 

This  last  formula  is  a  generalization  of  the  standard 
stationary  phase  formula  corresponding  to  H(s )  =  exp(±is). 
(See  Example  1  below.) 

To  illustrate  what  happens  when  logarithms  appear 
in  the  expansion  of  G(s),  as  s  -  0+,  let  us  suppose 
that  c  =  a  in  (1.2),  cp  is  as  in  (5*2),  and 


(5-13)  g(t)  ~  g01(t-c)iJ”1log(t-c)  +  g00(t-c)^1,  t  -  c+. 

After  some  calculation  we  find  that,  in  (5-15).  N( 0 )  =  1 
and 


(5-14) 


J00 


*01,  ,-”Ao 

V 

r0 

gpllogl  up  I  ~l 


Thus,  it  follows  from  (4.2)  and  (5.14)  that,  in  this  case. 
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As  a  final  general  result,  let  us  obtain  an  infinite 
asymptotic  expansion  of  Ic(X)  in  the  case  where 

(5.16)  Ci(t)  s  t-c 


and  g(t)  is  infinitely  differentiable  at  t  «*  c.  Then 


,(m) 


g(t)  ~  y  &-lri£l(t-c)m/  t  -  c+ 


(5.17) 


m*0 


g 


<t>  ~  y .(-hy t.,_. 


m=0 


Now  upon  applying  (4.11)  we  obtain 
®  ,-(m+l)  (m) ,  > 

(5.18)  Ic(>0  ~  2_  - - A£l[M[H(s)  ;m+l]  +  (-l)mMfH(-s ) ;m+li } . 

m=0 

If  any  terms  in  this  sum  (5-18)  are  non-zero,  then 
we  must  conclude  that  t  ■=  c  is  a  critical  point  for 
I(x).  Alternatively,  if  the  right  hand  side  of  (5.18) 
is  identically  zero,  then  t  -  c  is  not  critical.  The 
issue  depends  solely  on  the  kernel  and  at  that  only 
through  the  quantities 


(5.19)  M[H(s);m+l]  +  ( -l)mMrH( -s ) ;m+13 ,  m  -  0,1,2, - 
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Thus,  in  general,  t  =  c  is  critical  whenever  H  is 
such  that  at  least  one  of  the  quantities  (5.19)  is 
non-zero.  Furthermore,  it  is  readily  seen  that  the 
same  conclusion  holds  when  n  is  any  c“  function  that 
vanishes  at  t  =  c. 

We  shall  now  consider  two  illustrative  examples. 

Example  1:  Suppose  H(s)  is  the  complex  Fourier  kernel 
exp(is).  We  have  by  direct  computation 

(5.20)  M[exp(is ) ;z]  =  r(z)exp(2|^) 
and 

(5.21)  M[exp(is);z]  =  e“TTizM[exp(-is )  ;z] . 

From  this  last  relation  we  find 

(5.22)  M[exp(is ) ;m+l]  +  (-l)mM[exp(-is) ;m+l]  =0,  m  =  0,1,2,. 

and  hence,  as  anticipated  in  the  introduction,  the 
interior  zeros  of  are  not  critical  points  for  Fourier 
type  integrals. 

Suppose  now  that  g  is  continuous  at  t  >=  c  and  cp 
has  a  simple  stationary  point  there.  Then  it  follows 
from  (5.12)  and  (5.20)  that 

(5-23)  ICU)  ~  g(  c )  A^^-p4-^-|-cxp(sgn  rp"  ( c  )^p|. 

This  will  be  recognized  as  the  standard  stationary  phase 
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formula  In  the  case  where  ©(c)  -  0.  This  last  restric¬ 
tion  is  of  course  unnecessary  and  can  be  avoided  quite 
simply.  Indeed,  suppose  that  at  t  *  c,  ©  has  a 
simple  stationary  point,  but  ©(c)  /  0.  Then  we  write 

»b 

(5*24)  I(x)  =  exp(iX©(c))  exp[iX(©(t)-©(c)]g(t)dt. 

J  a 

Since 


(5.25) 


*  -  c(t)-©(c) 


has  a  sin.ple  stationary  point  at  t  *  c  and  $(c)  =  0, 
we  find  that  we  need  only  multiply  (5-23)  by  exp{iX©(c)} 
to  obtain  the  valid  result  in  this  case.  Furthermore, 
the  contribution  from  any  critical  point  at  which  ©  4  0 
can  be  recovered  in  an  analogous  manner  from  the  corre¬ 
sponding  contribution  in  the  case  where  ©  vanishes  at 
the  critical  point. 

Finally  suppose  that  c  =  a  and  (5-15)  holds.  Sup¬ 
pose  further  that  ©(t)-©(c)  satisfies  the  relations 
satisfied  by  ©(t)  in  (5-2).  Then  from  (5.15),  (5.20) 
and  the  remarks  of  the  preceeding  paragraph,  we  find 
that  now 


IC(M 


-f.V108’' +  T^l0Bl  "o'  -  Boo  -  +  Th)}- 


(5.26) 
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Here  f(z)  is  the  logarithmic  derivative  of  the  gamma 
function  F(z).  We  might  point  out  that  the  last  case 
was  considered  in  detail  by  Erdelyi  [  8 ]  and  by  r^Kenna  [  9 ] . 


Example  2\  Let  us  now  suppose  that  in  (1.2) 


(5.27)  H(s )  =  Jn(s),  n  >=  0,1,2, - 


Here  Jn  is  the  Bessel  function  of  the  first  kind  of 
order  n.  We  have  [ 5  ] 

_  1  r(i z  + 

(5-28)  M[  J  ( s ) ; z ]  =  2Z  x  rd-  - - . 

n  r(?n  “  ?z+i) 

Since  Jn(s)  is  even  about  zero  when  n  is  even  and 
odd  about  zero  when  n  is  odd,  we  have 

(5.29)  M[Jn(s);z]  *  -M(Jn(-s)jz],  n  odd 

M[Jn(s)jz]  *  M[Jn(-s);z),  n  even. 

Prom  this  it  follows  that 

(5.50)  M[Jn(s);m+l]  +  (-l)mM[Jn(-s)ja+l]  =  --  - M”*"* ■  - ' * [  1+  ( -1 ) m+] n 
n  n  r(^-[n+l-m] ) 

m  =  0,1,2, . 


and  hence  for  any  integer  n  one  half  of  the  quantities 
(5.30)  are  not  zero.  Thus  t  =  c  is  a  critical  point 


in  this  case. 
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A1 


APPENDIX 

In  this  Appendix,  we  shall  prove  Lemmas  1,  2  and  6 
of  the  text. 


Proof  of  Lemma  _1:  We  introduce  the  functions 

M  N(m)  _r 

(Al)  ck(s)  =  exp{-s‘k+isy-o(s)]-  £  V  Cjnns  m(log  s)n 

m=0  n=0 

( A2 )  Hk(s)  =  H(s )  -  ^k(s). 


Here 

(A3) 


u>(s) 


Ks>/  6 


I 


b 


e 


-6  J 


and  for  any  positive  k,  we  choose  M  =  M(k)  to  be 
the  largest  integer  such  that 


(A4)  Re(rM-rQ)  <  k. 

We  observe  that  H(s)  and  ak  have  identical 

asymptotic  expansions,  as  s  -  +®,  to  order 
"rM  -k+Re(r  ) 

s  M(log  s)"'  As  a  result,  Hk(s)  =  0(s  )  from 

which  it  follows  that  MfHk(s);zl  is  analytic  in  a 

strip  with  right  limit  Re(z)  <  k+Re(rQ).  The  real  exponential 

factor  in  ak  assures  us  that  Mr(7k(s);z-'  is  analytic 

in  the  left  half  plane  Re(z)  <  Re(rQ).  Thus  the  left 

limit  of  the  strip  of  analyticity  of  MrHk;z}  is  the 

same  as  that  of  MrH;z].  Let  us  denote  this  limit  by 

Re  ( z  )  =  a . 


A2 


Below,  we  list  the  relevant  Mellln  transforms  along 
with  their  strips  of  analyticity.  We  recall  that  in  its 
strip  of  analyticity  a  Mellin  transform  decays  to  zero 
as  ( y |  -  ®. 


Mellin  Transform  Strip  of  Analyticity 

M[H;z]  -a  <  Re(z)  <  Re(rQ) 

MCo^jz]  -®  <  Re(z)  <  Re(rQ) 

MrHR;z]  -a  <  Re(z)  <  Re(rQ)  +  k 

Prom  the  above  list  we  see  that  in  order  to  analytically 
continue  MfH;z]  into  the  region  Re(z)  <  k+Re(rQ)  we  need  only 
determine  the  analytic  continuation  of  M[ck;z]  into  the 
strip  Re(rQ)  <  Re(z)  <  Re(rQ;  +  k.  Furthermore,  since 
M[Hk?z]  decays  to  zero  as  |y|  -  ®  in  -a  <  Re(z )  <  Re(rQ)  +  k, 
any  algebraic  growth,  in  this  limit,  ’  M[H;z]  must 
arise  from  M[rk;z]. 

Since  ^(s)  iB  a  finite  sum,  so  is  its  Mellin 
transform.  A  typical  term  in  M[<?k;z]  is  given  by 

(A5)  I[z;r]  =  J  exp|-s_ic+is  u)(s)|(log  s )nsz“r-1ds ,  Re(z)<Re(r). 

In  (A5)  we  rotate  the  path  of  Integration  onto  the  ray 
arg  s  =  9  where  0  <  9sgn(b0)  <  n/2v.  The  effect  of 
this  rotation  is  to  introduce  sufficient  decay  at  ® 
so  that  the  integral  in  (A5)  converges  for  all  z. 

Hence  I[z;r]  can  be  continued  into  the  entire  z-plane 


as  a  holomorphic  function.  It  remains  only  to  estimate 


the  continuation  as  |yj  -  ®. 

In  (A5)  we  stretch  the  integration  variable  s  by 
the  factor  |  y  |  to  obtain 

(A6)  I[z;r]  =  |y|(x-r)/vj"exp{-(iy|1/vs)-k}sx-^1 

•  exp|i|y|  (svao(sly|1/v)  +  sgn  y  log  sjds. 

The  function  -u(  s  |  y |  is  a  finite  sum  of  the  form 

(A7)  '^(s|y|1>/v)  =  bQ  +  b1(sjy) 1/v)-6  +  .... 

From  this  we  find  that  the  phase  function 

( Afi  )  «(s,y)  =  sva)(sjy|  1/^v)  +  sgn  y  log  s 

has  stationary  points  (points  at  which  *g  =  0)  s^  such 
that 

<A9)  Sq  .  (^LZ)1/V+  0(|y|-6). 

Here  q  labels  the  different  choices  of  the  root. 

We  note  that  when  -sgn  y/vbQ  <  0,  no  stationary 
points  are  near  the  positive  real  axis.  When 
-sgn  y/vbn  "  0,  however,  there  are  simple  stationary 
points  on  or  near  the  positive  real  axis.  Armed  with 
this  information,  we  can  apply  the  results  of  reference  [7"'. 
to  conclude  that 


A4 


(A10) 

lim  i[z;rl  = 

|y|- 


0(  l y | 

o(  )y| 


"R),  for  all 

(x-r)/ v  -  \ 

2); 


R; 


zmu.  < 

vb0 


o 


o. 

vbo 


This  completes  the  proof. 

We  remark  that  the  relevant  results  of  reference  r  7  ] 
essentially  justify  the  formal  application  of  the  ordinary 
method  of  stationary  phase  to  the  integral  in  (A6). 


Proof  of  Lemma  2i  By  hypothesis 

(All)  M[ G ; z ]  «  f  siy"1(sxG(s))ds 

J0 

is  absolutely  convergent  for  all  Re(z)  =  x  >  xQ.*  Upon 

integrating  by  parts  q  times  and  using  the  stated 
/  d  \P  x 

properties  of  ( s )  (s  G) ,  we -obtain 

(A12)  M[G;z]  =  (^yj  jV^fs^ls^s )  )ds . 

We  now  break  the  interval  of  integration  at  the  points 
of  discontinuity  of  G^q+1^(s)  and  integrate  by  parts 
once  more  in  the  resulting  finite  sum  of  integrals,  in 
this  manner  we  obtain 

(A13)  f  siy'1(s~)q(sxG(s))ds  «  OdyT1) 

.  q 

Note  that  the  assumptions  made  imply  that  G(s)  -=  o(s  ), 

-xQ 

as  s  -  ®,  for  all  r  and  G{s)  =  0(s  )  as  s  -  0+. 


which,  when  combined  with  (A12),  completes  the  proof. 


Proof  of  Lemma  6:  If  Re(z)  <=  x  >  -Re(aQ),  then  the 
result  follows  from  Lemma  2  when  we  note  that  the  con¬ 
ditions  (sXG)  =  0,  p  *  0,1,..., q,  as  s  -  04- 

are  implied  by  the  assumed  differentiability  properties 
of  the  expansion  (2.15).  Now  suppose  that  o  is  any 
real  number  greater  than  Re(aQ).  Also  let  u(p)  and 
v(p)  be  positive  integers  satisfying  the  conditions 


(A14)  Re(a)i_1)  <  p  <  Re(a^),  Re(a0)  +  v  >  Re(a^). 


We  now  consider  the  functions 

A>-1  N(m)  a  ,x  v 

(A15)  -„(»)  -  (l  I  Vs  B(log  s)")e’B 

m=0  n=0 

(Al6)  Gp(s)  =  G(s )  -  ^p(s) 


and  note  that  v(o)  has  been  chosen  so  that  the  asymptotic 
expansions  of  G  and  agree  to  order  Re(au  Thus, 
as  s  -  0+ 


( A17 ) 


We  also  note  that  G  (s)  has  all  of  the  properties 

0 

attributed  to  G(s)  in  the  statement  of  the  lemma. 
Hence,  upon  applying  Lemma  2,  we  immediately  find  that. 


(Al8) 


M[Gp(s);zT  -  0(  |y|  • 


By  direct  calculation  we  have  that 
u-1  N(m)  m  a  +z-1 

(A19)  MtV*1  -  I  £d„Js"  (logs)V8 


m=0  n=0 

p.-l  N(m)  a  +z-l 

\  y  d  14  s  m 
L  L  mnrt  n\J 

m=0  n=0  az  0 


-sv 
e  ds 


y-1  N(m)d  +z  \ 

m_Q  n»0 


v 

ds 


in  the  region  Re(z)  >  -Re(aQ)  and  by  analytic  continu¬ 
ation  in  the  entire  z-plane.  We  know,  moreover,  that 
each  term  in  (A19)  decays  exponentially  as  |y|  -  ®  for 
all  x.  Finally,  since 


(A20)  M[G;z1  =  M(>  jz]  +  MfGp;z], 

we  have  that  (3*16)  holds  for  Re(z)  >  -Re^^  j).  How 
ever,  p  is  arbitrary  and  lim  Refa  ,  \1  =  «,  so  that 

p-«o  U  V  D  y 

upon  letting  p  -  ®  we  obtain  the  desired  result. 


